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Abstract 

We describe our current understanding on the phase transition phenomenon of the 
graph Laplacian eigenvectors constructed on a certain type of unweighted trees, 
which we previously observed through our numerical experiments. The eigen- 
value distribution for such a tree is a smooth bell-shaped curve starting from the 
eigenvalue up to 4. Then, at the eigenvalue 4, there is a sudden jump. Interest- 
ingly, the eigenvectors corresponding to the eigenvalues below 4 are semi-global 
oscillations (like Fourier modes) over the entire tree or one of the branches; on 
the other hand, those corresponding to the eigenvalues above 4 are much more lo- 
calized and concentrated (like wavelets) around junctions/branching vertices. For 
a special class of trees called starlike trees, we obtain a complete understanding 
of such phase transition phenomenon. For a general graph, we prove the number 
of the eigenvalues larger than 4 is bounded from above by the number of vertices 
whose degrees is strictly higher than 2. Moreover, we also prove that if a graph 
contains a branching path, then the magnitudes of the components of any eigen- 
vector corresponding to the eigenvalue greater than 4 decay exponentially from 
the branching vertex toward the leaf of that branch. 
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1. Introduction 



In our previous report [fTTI. we proposed a method to characterize dendrites of 
neurons, more specifically retinal ganglion cells (RGCs) of a mouse, and cluster 
them into different cell types using their morphological features, which are derived 
from the eigenvalues of the graph Laplacians when such dendrites are represented 
as graphs (in fact literally as "trees"). For the details on the data acquisition and 
the conversion of dendrites to graphs, see [ITTI and the references therein. While 
analyzing the eigenvalues and eigenvectors of those graph Laplacians, we ob- 
served a very peculiar phase transition phenomenon as shown in Figure [TTT| The 
eigenvalue distribution for each dendritic tree is a smooth bell-shaped curve start- 
ing from the eigenvalue up to 4. Then, at the eigenvalue 4, there is a sudden 
jump as shown in Figure [TTTJc, d). Interestingly, the eigenvectors corresponding 
to the eigenvalues below 4 are semi-global oscillations (like Fourier cosines/sines) 
over the entire dendrites or one of the dendrite arbors (or branches); on the other 
hand, those corresponding to the eigenvalues above 4 are much more localized 
and concentrated (like wavelets) around junctions/branching vertices, as shown 



in Figure 1.2 



We want to answer the following questions: 

Ql Why does such a phase transition phenomenon occur? 

Q2 What is the significance of the eigenvalue 4? 

Q3 Is there any tree that possesses an eigenvalue exactly equal to 4? 

Q4 What about more general graphs that possess eigenvalues exactly equal to 4? 

As for Ql and Q2, which are closely related, we have a complete answer for a 
specific and simple class of trees called starlike trees as described in Section [3| 
and a partial answer for more general trees and graphs such as those representing 
neuronal dendrites, which we discuss in Section]?} For Q3, we identify two classes 
of trees that have an eigenvalue exactly equal to 4, which is necessarily a simple 
eigenvalue, in Section [5j 

In Section [6| we will prove that the existence of a long path between two 
subgraphs implies that the eigenvalues of either of the subgraphs that are larger 
than 4 are actually very close to some eigenvalues of the whole graph. Then, 
in Section [7} we will give a counterexample to the conjecture that the largest 
component in the eigenvector corresponding to the largest eigenvalue (which is 
larger than 4) lies on the vertex of the highest degree. Finally, we describe our 
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Figure 1.1: Typical dendrites of Retinal Ganglion Cells (RGCs) of a mouse and the graph Lapla- 
cian eigenvalue distributions, (a) 2D projection of dendrites of an RGC of a mouse; (b) that of 
another RGC revealing different morphology; (c) the eigenvalue distribution of the RGC shown 
in (a); (d) that of the RGC shown in (b). Regardless of their morphological features, a phase 
transition occurs at the eigenvalue 4. 
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(a) RGC#100; Auai = 3.9994 



(b) RGC #100; /lii42 =4.3829 



Figure 1.2: The graph Laplacian eigenvectors of RGC #100. (a) The one corresponding to the 
eigenvalue An4i = 3.9994, immediately below the value 4; (b) the one corresponding to the 
eigenvalue /I1142 - 4.3829, immediately above the value 4. 

investigation on Q4 in Section [8] But let us first start by fixing our notation and 
reviewing the basics of graph Laplacians in Section |2] 



2. Definitions and Notation 

Let G = (V, E) be a graph where V = V{G) = {vi , V2, . . . , v„} is a set of vertices 
in G and E = E(G) = {61,62, ■■■ , 6,,,} is a set of edges where connects two 
vertices v,-, Vj for some 1 < i, j < n, and we write 6^ = (v,, Vj). Let du = divk) 
be the degree of the vertex Vk- If a graph G is a tr66, i.e., a connected graph 
without cycles, then it has m = n - \ edges. Let L(G) := D{G) - A(G) be the 
Laplacian matrix where D{G) := diag((ii, . . . ,dn) is called the d6gr66 matrix of 
G, i.e., the diagonal matrix of vertex degrees, and A(G) = (a,y) is the adjac6ncy 
matrix of G, i.e., aij = 1 if v, and Vj are adjacent; otherwise it is 0. Furthermore, 
let = Aq{G) < Xi{G) < ■■■ < A„-i{G) be the eigenvalues of L(G), and mciX) be 
the multiplicity of the eigenvalue A. More generally, if / c M is an interval of the 
real line, then we define moil) := #{Ak{G) e I}. 

At this point we would like to give a simple yet important example of a tree and 



its graph Laplacian: a path graph consisting of n vertices shown in Figure 2.3 



The graph Laplacian of such a path graph can be easily obtained and is instructive. 
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Figure 2.3: A path graph P„ provides a simple yet important example. 
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The eigenvectors of this matrix are nothing but the DCT Type II basis vectors used 
for the JPEG image compression standard; see e.g., lfT4]| . In fact, we have 



Ak = 4 sin 



nk 
2n 



(2.1) 
(2.2) 



for = 0, 1, . . . , n - 1, where = {(pi^k, ■ ■ ■ , (PhmY is the eigenvector correspond- 
ing to Ak- From these, it is clear that for any finite n 6 N, /Imax = < 4, and 
no localization/concentration occurs in the eigenvector (or any eigenvector), 
which is simply a global oscillation with the highest possible (i.e., the Nyquist) 

frequency, i.e., = {{-XT' sin(^ (; - ^)))^^ 



3. Analysis of Starlike Trees 

As one can imagine, analyzing this phase transition phenomenon for compli- 
cated dendritic trees turns out to be rather formidable. Hence, we start our analysis 
on a simpler class of trees called starlike trees. A starlike tree is a tree that has 
exactly one vertex of degree higher than 2. Examples are shown in Figure [34} 

We use the following notation. Let 5(ni, ^2, . . . , Hjt) be a starlike tree that has 
k{> 3) paths (i.e., branches) emanating from the central vertex vi. Let the rth 
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(a) 5(2,2,1,1,1,1) 



(b) 1, 1, 1, 1, 1, 1, 1) a.k.a. comet 



Figure 3.4: Typical examples of a starlike tree. 



branch have vertices excluding vi. Let > n2 > ■ ■ ■ > n^. Hence, the total 

k 

number of vertices is n = 1 + ^nt- 

i=\ 

Das proved the following results for a starhke tree S(n-[ .n^) in ||3l: 



'^max - ^n-\ < ^ + 1 + 



k-\ 



On the other hand, Grone and Merris proved the following lower bound for a 
general graph G with at least one edge: 

'Imax > max J(v^) + 1. (3.4) 

\<i<n 

Hence we have the following 

Corollary 3.1. A starlike tree has exactly one graph Laplacian eigenvalue greater 
than or equal to 4. The equality holds if and only if the starlike tree is Ki i, = 
5(1,1,1), which is also known as a claw. 



Proof. The first statement is easy to show. The lower bound in p.4[ ) is larger than 
or equal to 4 for any starlike tree since maxi<^<„ d(Vj) = d(vi) > 3. On the other 
hand, the second largest eigenvalue A„_2 is clearly strictly smaller than 4 due to 
([33]). 

To prove the second statement about the necessary condition on the equality 



(the sufficiency is easily verified), first note that by p.4[ ), d(vi) = 3 is a necessary 
condition for /l^ax = 4. Let di denote the highest degree of such a starlike tree, 
i.e., di = d(vi). Since we only consider starlike trees, the second highest degree 
d2 must be either 2 or 1. Now, we use the following 



Theorem 3.1 (Das 2004, Let G = {V,E) be a connected graph and di J2 
where d\ and dj are the highest and the second highest degree, respectively. Then, 
^max(G) = di + dj if and only ifG is a star graph. 

If /Imax(G) = 4 and di = 3, then using this theorem, we must have ^2 = 1- 
Hence, G must be a star graph with di = 3 and d2 = I, i.e., G = ^13. □ 

As for the concentration/localization of the eigenvector corresponding to 
the largest eigenvalue A^-i, we prove the following 

Theorem 3.2. Let 4>n-i = (^i.^-i^""" where (pj^n-i is the value of the 

eigenvector corresponding to the largest eigenvalue An-i at the vertex vj, j = 
I, . . . ,n. Then, the absolute value of this eigenvector at the central vertex vi can- 
not be exceeded by those at the other vertices, i.e., 

I<^i,«-il > j = 2,...,n. 

To prove this theorem, we use the following lemma, which is simply a corol- 
lary of Gerschgorin's theorem [fTSl Theorem 1.1]: 

Lemma 3.1. Let A be a square matrix of size n x n, Ak(A) be any eigenvalue 
of A, and <pf^ = {(pi^k, ■ ■ ■ ,(pn,kf the corresponding eigenvector Let k* denote 
the index of the largest eigenvector component in <f>j^, i.e., \(f>k',k\ = maxygAr 
where N := {!,..., n}. Then, we must have Ak(A) e F^'iA), where FiiA) : = 
|z 6 C : \z- aul < YjjeN\{i\ is the ith Gerschgorin disk of A. In other words, 
for the index of the largest eigenvector component, the corresponding Gerschgorin 
disk must contain the eigenvalue. 



Proof. Recall the proof of Gerschgorin's theorem. The ^*th row of = A^ipi, 
yields 

jeN\{k'} \4'k*,k\ jeN\{k'} 



Zi I ki'-^i I I 
\ak'j\ i 1 < 2j r'^'jl ■ 



This implies Ak 6 A*(^), which proves the lemma. □ 



Proof of Theorem [372] First of all, by Corollary |3.1| we have A„-i > 4. However, 
An-i = 4 happens only for ^13. In that case, it is easy to see that this theorem holds 
by directly examining the eigenvector = ^3 oc (3,-1,-1,-1)^. Hence, let us 



examine the case > 4. In this case. Lemma 3.1 indicates 4 < A„-i e r(„_i).(L) 
where (n - I)* e N is the index of the largest component in Now, note that 
the disk r,(L) for any vertex v,- that has degree 2 is {z 6 C : |z - 2| < 2} (and 
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{z e C : |z - 1| < 1} for a degree 1 vertex). This means that the Gerschgorin 
disk r^n-iy containing the eigenvalue /l„_i > 4 cannot be in the union of the 
Gerschgorin disks corresponding to the vertices whose degrees are 2 or lower. 
Hence the index of the largest eigenvector component in must correspond to 
an index for which the vertex has degree 3 or higher. In our starlike-tree case, 
there is only one such vertex, vi, i.e., (n - 1)* = 1. □ 

For different proofs without using Gerschgorin's theorem, see Das [|3l Lemma 4.2] 
and E. Woei's dissertation [fT6l . We note that our proof using Gerschgorin's disks 
is more powerful than those other proofs and can be used for more general situa- 
tions than the starlike trees as we will see in Section |4l 

Remark 3.1. Let <p = {(pi,(p2, . . . , (pnf be an eigenvector of a starlike tree S (ni, . . . ,ni,) 
corresponding to the Laplacian eigenvalue A. Without loss of generality, letvj, ■ ■ . , v„j+i 
be the n^ vertices along a branch emanating from the central vertex Vi with v„,+i 
being the leaf ( or pendant) vertex. Then, along this branch, the eigenvector com- 
ponents satisfy the following equations: 

Mn+i = 0»i+i-^«p (3.5) 
A^j = 2(pj- (pj^i- (pj+i 2<j<ni. (3.6) 



From Eq. p.6[ ), we have the following recursion relation: 

(f)j+i +{A- 2)(pj = 0, J = 2, . . . , Hi . 

This recursion can be explicitly solved using the roots of the characteristic equa- 
tion 

r^ + {A-2)r+l= 0, (3.7) 



and when (3.7 1 has distinct roots r\,r2, the general solution can be written as 



(Pj = Ar{~^ + Br{~\ j = 2, . . . , m + 1, (3.8) 



where A, B are appropriate constants derived from the boundary condition p.5[ ). 
Now, let us consider these roots of p.7| ) in detail. The discriminant of p.7| ) is 



D{A) := (A-2f -4 = A{A-4). 

Since we know that A > 0, this discriminant changes its sign depending on A < 4 
or A > A. (Note that A = A occurs only for the claw Ki ^ on which we explicitly 
know everything; hence we will not discuss this case further in this remark.) If 
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A < 4, then D{A) < and it is easy to show that the roots are complex valued with 
magnitude 1. This implies that p.8| ) becomes 



(pj = A' cos(a>( J - 2)) + B' sinMj - 2)), j = 2, . . . , + 1 , 

where o) satisfies tan oj = V^^(4^^^/(2 - A), and A', B' are appropriate constants. 
In other words, if A < 4, the eigenvector along this branch is of oscillatory na- 
ture. On the other hand, if A > 4, then !D{A) > and it is easy to show that 
both ri and r2 are real valued with -1 <ri = ^2 - A + y/A(A - 4)j /2 < while 

r2 = (2 - A - y/A(A — 4)j /2 < -1. On the surface, the term Br!^^ looks like a 
dominating part in p.SP ; however, we see from p.5| ) that 10,,, | > |0„j+i|, which 



means the real dominating part in p.8| ) for j = 2, . . . ,n\ + \ is the term Ar^^'^. 
Hence we conclude that \(pj\ decays exponentially with j, that is, the eigenvector 
component decays rapidly towards the leaves. The siuation is the same for the 
other branches. 

In summary, we have shown that a starhke tree has only one eigenvalue > 4, 
and its eigenvector is localized at the central vertex in the sense of Theorem 3.1. 
Furthermore, the other eigenvectors are of oscillatory nature. Therefore the phase 
transition phenomenon for a starlike tree (with the eigenvalue 4 as its threshold) 
is completely understood. 

4. The Localization Phenomena on General Graphs 

Unfortunately, actual dendritic trees are not exactly starlike. However, our 
numerical computations and data analysis on totally 179 RGCs indicate that: 

#{ J 6 |J(v,) > 2} - mG([4, CX3)) 

< — < 0.047 

n 

for each RGC. Hence, we can define the starlikeliness S£{T) of a given tree T as 

^^^^^ ^ _ #{jeA^|^(vy)>2}-mr([4,cx,)) 

n 

We note that S£{T) = 1 for a certain class of RGCs whose dendrites are sparsely 
spread (see [11] for the characterization). This means that dendrites in that class 
are all close to a starlike tree or a concatenation of several starlike trees. We 
show some examples of dendritic trees with S£{T) = 1 and with S€{T) < 1 in 



Figure 4.5 



The above observation has led us to prove the following 
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(a) RGC#100;5^(r)= 1 



(b) RGC #155; S{{T) = 0.953 < 1 



Figure 4.5: Zoomed-up versions of parts of some dendritic trees. 



Theorem 4.1. For any graph G of finite volume, i.e., 2"=i d{vj) < oo, we have 

0<mG([4,cx,))<#{jeA^|J(vy)>2} 

and each eigenvector corresponding to A> 4 has its largest component ( in abso- 
lute value ) on the vertex whose degree is higher than 2. 

We refer the interested readers to [fTOl Sec. 2] that reviews various relationships 
between the multiplicity of certain eigenvalues and the graph structural properties 
different from our Theorem 14. 1[ 



Proof. The second statement follows from Lemma 3.1 because the Gerschgorin 



disks corresponding to vertices of degree 1 or 2 do not include A > A. 

We next prove the first statement. Let L be a Laplacian matrix of G. We can 
apply a permutation P such that 



P^LP 



E L2 



(4.9) 



where the diagonals of Li are 3 or larger (correspond to vertices of degree > 2), 
and the diagonals of L2 are 2 or 1. Suppose L2 is ^-by-^. By Gerschgorin's theorem 
all the eigenvalues of L2 must be 4 or below. 

In fact, we can prove the eigenvalues of L2 are strictly below 4. By ifTSl The- 
orem 1.12], for an irreducible matrix (a Laplacian of a connected graph is irre- 
ducible) an eigenvalue can exist on the boundary of the union of the Gerschgorin 
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£ >3 (for example when £ = 3, L2 



disks only if it is the boundary of all the disks. Furthermore, if there is such an 
eigenvalue, then the corresponding eigenvector has the property that all its com- 
ponents have the same absolute value. 

Suppose on the contrary that L2X = Ax. Suppose without loss of generality 
that L2 is irreducible; if not, we can apply a permutation so that PL2P^ is block 
diagonal and treat each block separately. 

Now if L2 has a diagonal 1 , then the corresponding Gerschgorin disk lies on 
[0, 2], which does not pass 4. Hence by [15, Theorem 1.12] this case is ruled out. 
It follows that all the diagonals of Li are 2, and the sum of the absolute values of 

the rows of L2 are all 4 (this happens only if L2 is disjoint from Li). So we need 

2 -1 -11 

-1 2 -1 ). Now the zth (1 < z < i) row 

-1-1 2_ 

of L2X = Ax and the fact \xi\ = \x2\ = ■ ■ ■ = \x[\ force x,- = -xj for all j 4^ i. 
This needs to hold for all z, which clearly cannot happen for ^ > 3. Therefore the 
eigenvalues of L2 must be strictly below 4. 

By the min-max characterization of the eigenvalues of P^LP, denoting by 
A({P^LP) the ^th smallest eigenvalue, we have 

AfiP'^LP) = min max v"^(P^LP)v. 

dim5=f>'espaii(5'),|l3'||2=l 

Hence letting 5 be the last € column vectors of the identity /„ and noting S ^P^LPS 
L2, we have 

Ae(P''LP)< max y\P''LP)y = A^.^S^P'^LPS,) 

3'eSo,lbll2=l 

= ^max('^'2)- 

Since /^max(^2) < A, we conclude that P^LP (and hence L) has at least £ eigen- 
values smaller than 4, i.e., mG([0,4)) > L Hence, mcdA, 00)) = n - mciiOA)) < 
n - £ = #{j e N\ d{vj) > 2}, which proves the first statement. □ 

To give a further explanation for the eigenvector localization behavior ob- 
served in Introduction, we next show that eigenvector components of A > A must 
decay exponentially along a branching path. 

Theorem 4.2. Suppose that a graph G has a branch that consists of a path of 
length k, whose indices are {z'l, Z2, . . . , ik) where z'l is connected to the rest of the 
graph and i^ is the leaf of that branch. Then for any eigenvalue A greater than 4, 
the corresponding eigenvector <f) = (0i, • • • ,0,,)^ satisfies 



ij+i 



<7\cpi,\ forj=\,2,...,k-l, (4.10) 



11 



where 



r := 



A-2 



< 1. 



(4.11) 



Hence \(pi\ < y-'~^\(f>i^\for j = 1,. . . ,k, that is, the magnitude of the components of 
an eigenvector corresponding to any A > 4 along such a branch decays exponen- 
tially toward its leaf with rate at least y. 

Proof. There exists a permutation P such that 

L, £T 



L := P^LP 



E U 



where 



2 -1 
-1 2 -1 

-1 2 -1 



ikxk 



-1 2 -1 

-1 1 _ 

and E has a -1 in the top-right comer and elsewhere. The diagonals of L2 
correspond to the vertices v,p . . . , v,^ of the branch under consideration. 

Let = A<l) with A> A. We have Ly = Ay where J = (3^1,^2, •• • , Jn)^ = P^^- 
Note that Cy„_^+ 1 , yn-t+i , • • • , j„ ) = (0,-, , , . . . , 0,-, ) . The last row of Ly = Ay gives 



hence 



\yn\ 



-yn-\ +yn = Ayn, 

-\yn-i\ < ybn-il 



A- 1 

The {n - l)st row of Ly = Ay gives 

-yn-2 + 2y„-i -yn = Ay^-i. 

Using \yn\ < \y„.i\ we get 

I I _ \yn-2+yn\ ^ \yn-2\ + bn-ll 

- - A-2 ' 

from which we get \yn^i\ < \yn-2\- Therefore ly^l < \y„-i\ < \y„-2\, and so 

\yn-2 + yn\ ^ 'Ayn-i\ 



(4.12) 



(4.13) 



\yn- 



A-2 



< 



A-2 



y\yn- 



21- 



Repeating this argument k-\ times we obtain (4.10). 



□ 
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We note that the inequalities ( |4.12[ ) and ( |4.13[ ) include considerable overes- 
timates, and tighter bounds can be obtained at the cost of simplicity. Hence in 
practice the decay rate is much smaller than y defined in ( |4.1 1[ ). We also note that 
the larger the eigenvalue A > 4, the smaller the decay rate y is, i.e., the faster the 
amplitude decays along the branching path. 

Also note that the above result holds for any branching path of a tree. In partic- 
ular, if a tree has k branches consisting of paths, they must all have the exponential 
decay in eigenvector components if A > 4. This gives a partial explanation for the 
eigenvector localization behavior observed in Introduction. However, the theorem 
cannot compare the eigenvector components corresponding to branches emanat- 
ing from different vertices of degrees higher than 2, so a complete explanation 
remains an open problem. 

Remark 4.1. Let us briefly consider the case A = 4. In this case we have 7 = ^ = 
1, suggesting the corresponding eigenvector components along a branching path 
may not decay. However, we can still prove that unless (pi^ = (pi^_ = ■ ■ ■ = (pii^ = 0, 
we must have 

|04l<|0„_J<---<|0,-J. (4.14) 

In other words, the eigenvector components must decay along the branch, al- 
though not necessarily exponentially. To see this, weflrst note that ify^ = 0, then 
the last row of Ly = Ay forces y„^i = 0. Then, j„ = = together with the 
(n - \)st row gives yn-i - 0. Repeating this argument we conclude that yj must be 
zero for all j = n - k + 1, . . . ,n. Now suppose that > 0. Following the above 
arguments we see that the inequality in ( |4.12[ ) with y = I must be strict, that is, 
\yn\ < \yn- il Using this we see that the inequality in ( |4.13[ ) must also be strict, 
hence |y„_i| < \yn-2l Repeating this argument proves (|4.14|). 



5. A Class of Trees Having the Eigenvalue 4 

As raised in Introduction, we are interested in answering Q3: Is there any tree 
that possesses an eigenvalue exactly equal to 4? To answer this question, we use 
the following result of Guo [9] (written in our own notation). 

Theorem 5.1 (Guo 2006, ['91). Let T be a tree with n vertices. Then, 



Aj(T) < 



j = 0,...,n- I, 



and the equality holds if and only if all of the following hold: a) j Q; b) n — j 
divides n; and c) T is spanned by n - j vertex disjoint copies of K, 

' I'-j 
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(a) 



(b) 



Figure 5.6: (a) A tree spanned by multiple copies of K[ y, connected via their central vertices. This 
tree has an eigenvalue equal to 4 with multiplicity 1 . (b) The eigenvalue distribution of such a tree 
spanned by 5 copies of Ki y,. We note that StiJ) - 1 for this tree. 



Here, a tree T = T{V,E) is said to be spanned by i vertex disjoint copies of 
identical graphs Ki{Vi, £,) for / = !,...,£ if V = Uf=i Vi and V,- n = for all 
i 4^ j. Figure [53) ^a) shows an example of such vertex disjoint copies for Ki = Ki j, 
by connecting their central vertices. We note that there are many other ways to 
form disjoint vertex copies of Ki ^,. 
This theorem implies the following 

Corollary 5.1. A tree has an eigenvalue exactly equal to 4 if it is spanned by 
m(= n/4 e N) vertex disjoint copies ofKi j = 5(1, 1, 1). 



Figure 5.6 'b) shows the eigenvalue distribution of a tree spanned by m = 5 
copies of as shown in Figure [5^ a). Regardless of m, the eigenvector corre- 
sponding to the eigenvalue 4 has only two values: one constant value at the central 
vertices, and the other constant value of the opposite sign at the leaves, as shown 



in Figure 5.7 a). By contrast, the eigenvector corresponding to the largest eigen- 



value is again concentrated around the central vertex as shown in Figure 5.7 ^b) 



Theorem 5.1 asserts that a general tree T with n vertices can have at most 
ln/4\ Laplacian eigenvalues > 4. We also know by Theorem 2.1 of [8] that any 
tree T possessing the eigenvalue 4 must have multiplicity mj{A) = 1 and n = 4m 



for some m e N. Hence, Theorem 5]T]also asserts that trees spanned by m vertex 



disjoint copies of form the only class of trees for which 4 is the 3n/4(= 3m)th 
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(a) <f>i5 (b) 4>ig 



Figure 5.7: (a) The eigenvector corresponding to Ai^ = 4 in the 3D perspective view, (b) 
The eigenvector corresponding to the maximum eigenvalue /I19 = 7.1091, which concentrates 
around the central vertex. 



eigenvalue of T; in other words, trees in this class are the only ones that have 
exactly n/4(= m) eigenvalues > 4. 

Trees spanned by vertex disjoint copies of ^13, however, are not the only ones 
that have an eigenvalue exactly equal to 4. For instance. Example 2.9 of [|8l, 
which has n = 36 = 4 • 9 vertices and is called Z4 as shown in Figure |5.8[ is 



non-isomorphic to any tree spanned by 9 vertex disjoint copies of Ki^\ yet it has 
/I30 = 4 (but = 1 4). 

On the other hand, we have the following 

Proposition 5.1. Ifn < 11, any tree possessing an eigenvalue exactly equal to 4 
must be spanned by vertex disjoint copies ofKi j. 

Proof. First of all, 4 must divide n, hence n = Am with m = 1 or m = 2. If m = 1, 
then we know from Theorem |5TT] that Amax = /I3 < 4 and /I2 < 2. Hence, /I3 = 4 is 
the only possibility, and consequently T = K13 using the same theorem. If m = 2, 
then Theorem |5 . 1 1 states that A^^x = Aj < S; A(, < 4; and A5 < 3; . . .If A^ = 4, then 
the necessary and sufficient conditions for the equality in Theorem |5.1| state that 
T must be spanned by two vertex disjoint copies of ^13, so we are done. Now we 
still need to show that Aj cannot be 4. Let di be the degree of the highest degree 
vertex of a tree T under consideration. Then, we have 

di + \ < At < di +2^Jdl - I. (5.15) 

Here, the lower bound is due to Grone and Merris Q and the upper bound is due 
to Stevanovic [fT3l . Now, we need to check a few cases of the values of di. 
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Figure 5.8: Yet another tree Z4 (Example 2.9 of IHl with k - 4) that has an eigenvalue exactly 
equal to 4. This is non-isomorphic to any tree spanned by vertex disjoint copies of K12 such as the 



one shown in Figure 5.6 ^a). 



If di = 2, then the upper bound in ( |5.15[ ) is 4. Hence, Aj = 4 cannot happen. 
(This includes the case of a path graph that cannot reach the eigenvalue 4). 
So, we must have di >3. 



• If di > 3, of course, the lower bound in ( |5.15[ ) is greater than 4. Hence, Aj 
cannot be 4 either. 

• Finally, if di = 3, then the above bounds are: 4 < < 5.8284- • • . Can 
Aj = 4 m this case? According to Zhang and Luo [TTj, the equality in that 
lower bound holds if and only if there exists a vertex that is adjacent to all 
the other vertices in T. That is, the degree of that vertex is n - 1 = 7. Since 
di = 3, this cannot happen. 

Hence, for « = 8, the only possibility for a tree T to have an eigenvalue exactly 
equal to 4 is the case when A^ = 4, which happens if and only if T is spanned by 
two vertex disjoint copies of ^1 3. □ 

It turns out, however, that proving the necessity for n > 1 1 using similar argu- 
ments quickly becomes cumbersome, even for the next step n = 12. At this point, 
we do not know whether there are other classes of trees than Z4 discussed above or 
those spanned by vertex disjoint copies of ^13 that can have an eigenvalue exactly 
equal to 4. Hence, identifying every possible tree that has an eigenvalue exactly 
equal to 4 is an open problem. 
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6. Implication of a Long Path on Eigenvalues 



In Section|4]we saw that for a graph that has a branch consisting of a long path, 
its Laplacian eigenvalue greater than 4 has the property that the corresponding 
eigenvector components along the branch must decay exponentially. 

Here we discuss a consequence of such a structure in terms of the eigenvalues. 
We consider a graph G formed by connecting two graphs Gi and G3 with a path 
G2. Note that this is a more general graph than in Section|4] (which can be regarded 
as the case without G3). We show that if G2 is a long path then any eigenvalue 
greater than 4 of the Laplacian of either of the two subgraphs Gi U G2 and G2 U G3 
must be nearly the same as an eigenvalue of the Laplacian of the whole graph G. 

Theorem 6.1. Let G be a graph obtained by connecting two graphs with a path, 
whose Laplacian L can be expressed as 



L = 



L, 
El 




El 





£■2 



where E\ and E2 have -1 in the top-right corner and elsewhere. Li is ti x tifor 
i = 1,2,3 and L^ represents the path G2, that is, a tridiagonal matrix with 2 on 
the diagonals and -1 on the off-diagonals. 

Let A> Abe any eigenvalue of the top-left (£1 + £2) x (^1 + ^2) (or bottom-right 
{€2 + £3) X (£2 + £3)) submatrix of L. Then there exists an eigenvalue AofL such 
that 

\A-l\<y'\ (6.16) 

2 



where y := < 1. 

/i— 2 

Proof. We treat the case where X is an eigenvalue of the top-left {t\ + ^2) x < 
part of L, which we den ote by L12. The other case is analogous. 



+ ^2) 



As in Theorem 



4.2 



we can show that any eigenvalue X > A oi Lu has its 
corresponding eigenvector components decay exponentially along the path G2. 
This means that the bottom eigenvector component is smaller than y^^ in absolute 
value (we normalize the eigenvector so that it has unit norm) where 7 ■= < 1 
as in ( |4.11[ ). 

Let L\2 = QAQ^ be an eigendecomposition where Q = / and the eigen- 
values are arranged so that X appears in the top diagonal of A. For notational 
convenience let Cn '■= £\ + ^2- Then, consider the matrix 



(6.17) 








L 


Q 







A ve\ 





/ 







/ 




ey L3 
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Figure 7.9: Counterexample graph for the conjecture. 



where Ci = (1,0, ... ,0)^ 6 M^' ^nd v = (vi, . . .,ve^^y e M^'^ Direct calculations 
show that Vi = -qti2,i where qe^^^i is the bottom component of the eigenvector of 
L\2 corresponding to the /th eigenvalue. In particular, by the above argument we 

have|^f,,,il = Ivil <r^H« !)• 

Note that in the first row and column of L, the only nonzero s are the diagonal 
(which is A), and the (1, ^12 + 1) and {in + 1,1) entries, both of which are equal to 
vi. Now, viewing the (1, ^12 + 1) and {£12 + 1,1) entries of L as perturbations (write 
L = Li + L2 where Li is obtained by setting the (1, ^12 + 1) and (/'12 + 1, 1) entries 
of L to 0) and using Weyl's theorem [6, Theorem 8.1.5] we see that there exists an 
eigenvalue Aof L (and hence of L) that lies in the interval [A - IIL2II2, A + IIL2II2] = 
[A- |vi|,l+ |vi|]. Together with |vi| < we obtain (6.161). □ 



Recall that y'^^ decays exponentially with £2, and it can be negligibly small for 
moderate £2', for example, for (/I, ^2) = (5,30) we have y^- = 5.2 x 10"^. We 
conclude that the existence of a subgraph consisting of a long path implies that 
the eigenvalues /I > 4 of a subgraph must match those of the whole graph to high 
accuracy. 



7. On the Eigenvector of the Largest Eigenvalue 

In view of the results in Section |4] it is natural to ask whether it is always 
true that the largest component of the eigenvector corresponding to the largest 
eigenvalue of a Laplacian matrix of a graph lies on the vertex of the highest degree. 
Here we show by a counterexample that this is not necessarily true. 

Consider for example a tree as in Figure |7.9[ which is generated as follows: 
first we connect m copies of K12 (equal to ^3) as shown in Figure [5^a); then add 



to the right a comet 5(/', 1, 1, 1, 1) as in Figure 3.4 b). 

Now for sufficiently large m and £ (m, £ > 5 is sufficient), the largest compo- 
nent in the eigenvector ^ corresponding to the largest eigenvalue of the resulting 
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Laplacian L occurs at one of the central vertices of K12, not at the vertex of degree 
5 belonging to the comet. 

Let us explain how we came up with this counterexample. The idea is based 
on two facts. The first is the discussion in Section [6} where we noted that a long 
path G2 implies any eigenvalue larger than 4 must be close to an eigenvalue of 
a subgraph Gi U G2 or G2 U G3. Therefore, in the notation of Section |6} by 
connecting two graphs (Gi = mKi 2 and G3 = ^1,5, a star) with a path G2 such that 
the largest eigenvalue A of L12 is larger than that of L23, we ensure that the largest 
eigenvalue A of Lis very close to A. The second is the Davis-Kahan sin 6 theorem 
O, which states that a small perturbation of size the matrix Li (recall the 
proof of Theorem 6.1[ ) can only induce small perturbation also in the eigenvector: 



its angular perturbation is bounded by y^^/S, where 6 is the distance between A 
and the eigenvalues of L after removing its first row and column. Furthermore, 
the eigenpair (A, ^) of Li satisfies ^ = (1,0,..., 0)^, and the eigenvectors ^ (^ ^ 
by Davis-Kahan) of L and <p of L corresponding to A are related by ^ = [ q / ] ^' 



which follows from ( |6.17[ ). Therefore, ^ has its large components at the vertices 
belonging to G\. In view of these our approach was to find two graphs G\ and G3 
such that the highest degrees of the vertices of Gi and G3 are 4 and 5, respectively, 
and the largest eigenvalue of the Laplacian of Gi is larger than that of G3. 

8. Discussion 

In this paper, we obtained precise understanding of the phase transition phe- 
nomenon of the combinatorial graph Laplacian eigenvalues and eigenvectors for 
starlike trees. For a more complicated class of graphs including those representing 
dendritic trees of RGCs, we proved in Theorem |4.1| that the number of the eigen- 
values greater than or equal to 4 is bounded from above by the number of vertices 



whose degrees are strictly higher than 2. In Theorem 4.2, we proved that if a 
graph has a branching path, the magnitude of the components of an eigenvector 
corresponding to any eigenvalue greater than 4 along such a branching path de- 



cays exponentially toward its leaf. In Remark 4.1[ we also extended Theorem 4.2 



for the case of A = A although the decay may not be exponential. 

As for Q3 raised in Introduction — "Is there any tree that possesses an eigen- 
value exactly equal to 4?" — we showed that any tree with n vertices {n = Am for 
some m e N) spanned by m vertex disjoint copies of ^13 possesses an eigenvalue 



exactly equal to 4 in Corollary 5.1[ and that such class of trees are the only ones 



that can have an eigenvalue exactly equal to 4 if n < 11 in Proposition 5.1 On 



the other hand, for larger n, we pointed out that not only those spanned by vertex 
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disjoint copies of ^13, but also a tree called Z4 discovered in ^ and shown in 



Figure 5.8 have an eigenvalue exactly equal to 4. A challenging yet interesting 
question is whether or not one can identify every possible tree that has an eigen- 
value 4. 

Another quite interesting question is Q4 raised in Introduction: "Can a simple 
and connected graph, not necessarily a tree, have eigenvalues equal to 4?" The 
answer is a clear "Yes." For example, the J-cube (d > 1), i.e., the J-fold Carte- 
sian product of K2 with itself is known to have the Laplacian eigenvalue 4 with 
multiplicity d(d - l)/2; see e.g., [1, Sec. 4.3.1]. 

Another interesting example is a regular finite lattice graph in M'', d > I, which 



is simply the d-iold Cartesian product of a path shown in Figure 2.3 with itself. 
Such a lattice graph has repeated eigenvalue 4. In fact, each eigenvalue and the 
corresponding eigenvector of such a lattice graph can be written as 



d 



2nJ 



i=l 

I i .ttI v. 4- ^ 



(8.19) 



where j,, e Z/nZ for each /, as shown by Burden and Hedstrom ||2l. Note that 
( |8.18[ ) and ( |8.19[ ) are also valid for d = 1 . In that case these reduce to ( |2.2[ ) that 



we already examined in Section |2j 

Now, determining mdA-), i.e., the multiplicity of the eigenvalue 4 of this lattice 
graph, is equivalent to finding the number of the integer solutions (71, ... , j^) e 
(Z/nZY to the following equation: 

tsi„^(g)=l. (8.20, 

(=1 

For d = I, there is no solution as we mentioned in Section|2] For J = 2, it is easy to 
show that mo (4) = r - 1 by direct examination of ( |8.20[ ) using some trigonometric 



identities. For d = 3, mc(4-) behaves in a much more complicated manner, which 
is deeply related to number theory. We expect that more complicated situations 
occur for d > 3. We are currently investigating this on regular finite lattices. On 



the other hand, it is clear from ( |8.19[ ) that the eigenvectors corresponding to the 
eigenvalues greater than or equal to 4 on such lattice graphs cannot be localized or 
concentrated on those vertices whose degree is higher than 2 unlike the tree case. 
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Theorem 4.2 and Remark 4[T]do not apply either since such a finite lattice graph 
do not have branching paths. 

Finally, we would like to note that even a simple path, such as the one shown 



in Figure [23| exhibits the eigenfunction localization phenomena if it has nonuni- 
form edge weights, which we recently observed numerically. We will report our 
progress on investigation of localization phenomena on such weighted graphs at a 
later date. 



Acknowledgments 

We thank the referees for their remarks and suggestions. This research was 
partially supported by the following grants from the Office of Naval Research: 
NOOO 14-09- 1-0041; N00014-09-1-0318. A preliminary version of a part of the 
material in this paper [12] was presented at the workshop on "Recent development 
and scientific applications in wavelet analysis" held at the Research Institute for 
Mathematical Sciences (RIMS), Kyoto University, Japan, in October 2010, and at 
the 7th International Congress on Industrial and Applied Mathematics (ICIAM), 
held in Vancouver, Canada, in July 201 1. 



References 

[1] Biyikoglu T., Leydold J., Stadler R F. (2007) Laplacian Eigenvectors of 
Graphs. Lecture Notes in Mathematics 1915. Springer, New York. 

[2] Burden R. L., Hedstrom G. W. (1972) The distribution of the eigenvalues of 
the discrete Laplacian. BIT 12:475-488. 

[3] Das K. C. (2007) Some spectral properties of the Laplacian matrix of starlike 
trees. Italian J. Pure Appl. Math. 21:197-210. 

[4] Das K. C. (2004) A characterization on graphs which achieve the upper 
bound for the largest Laplacian eigenvalue of graphs. Linear Algebra Appl. 
376:173-186. 

[5] Davis C, Kahan W. M. (1970) The rotation of eigenvectors by a perturba- 
tion. III. SIAMJ. Numer Anal. 7:1-46. 

[6] Golub G. H., Van Loan C. F. (1996) Matrix Computations. The Johns 
Hopkins Univ. Press, Baltimore, MD, 3rd edition. 



21 



[7] Grone R., Merris R. (1994) The Laplacian spectrum of a graph II. SI AM J. 
Discrete Math. 7:221-229. 

[8] Grone R., Merris R., Sunder, V. S. (1990) The Laplacian spectrum of a 
graph. SIAM J. Matrix Anal. Appl. 11:218-238. 

[9] Guo J. M. (2007) The Ath Laplacian eigenvalue of a tree. / Graph Theor 
54:51-57. 

[10] Merris, R. (1994) Laplacian matrices of graphs: A survey. Linear Algebra 
Appl. 197/198:143-176. 

[1 1] Saito N., Woei E. (2009) Analysis of neuronal dendrite patterns using eigen- 
values of graph Laplacians. J SIAM Letters 1:13-16. Invited paper. 

[12] Saito N., Woei E. (2011) On the phase transition phenomenon of graph 
Laplacian eigenfunctions on trees. RIMS Kokyuroku 1743:77-90. 

[13] Stevanovic, D. (2003) Bounding the largest eigenvalue of trees in terms of 
the largest vertex degree. Linear Algebra Appl. 360:35^2. 

[14] Strang G. (1999) The discrete cosine transform. SIAM Review 41:135-147 . 

[15] Varga R. S. (2004) Gersgorin and His Circles, Springer, New York. 

[16] Woei E. (2012) Ph.D. dissertation. Dept. Math., Univ. CaUfomia, Davis. In 
preparation. 

[17] Zhang, X.-D., Luo, R. (2002) The spectral radius of triangle-free graphs. 
Australas. J. Combin. 26:33-39. 



22 



